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Abstract 

We study transport through a ferromagnetic single-electron transistor. The 
resistance is represented as a path integral, so that systems where the tun- 
nel resistances are smaller than the quantum resistance can be investigated. 
Beyond the low order sequential tunneling and co-tunneling regimes, a large 
magnetoresistance ratio at sufficiently low temperatures is found. In the op- 
posite limit, when the thermal energy is larger than the charging energy, the 
magnetoresistance ratio is only slightly enhanced. 
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Electron tunneling between two ferromagnetic metallic leads separated by an insulating 
layer has been studied since the seventies [IT]. The novel feature of these structures is that 
the conductance depends on the relative orientation of the magnetizations in the leads. 
By tuning the external magnetic field, the magnetization of a ferromagnetic metal can be 
controlled, and hence the conductance of the device changes. The influence of the applied 
magnetic field on the transport properties is characterized by the tunnel magnetoresistance 
ratio (TMR): 

7 = {Rap - Rp)/Rp , (1) 

where Rp {Rap) is the resistance when the magnetizations in the leads are parallel (anti- 
parallel). Devices with large TMR can be used in future magneto- resistive sensors or mag- 
netic random access memories [^. 

In the past years, reproducible measurements of TMR have been reported in struc- 
tures with ultrasmall tunnel junctions. When the capacitances of the tunnel junctions are 
so small that the charging energy is larger than the thermal energy, the current flow is 
blocked when the bias voltage is smaller than a critical value determined by the junction 
capacitances. Such a Coulomb blockade of electron tunneling through insulating barriers 
has been widely investigated when the leads are nonmagnetic ||7|-[T0|. The single-electron 



transistor consists of a metallic island connected to two metallic leads via tunnel junctions 
and coupled to a gate- voltage Ug via a capacitor Cg. The lowest order perturbation in terms 
of the dimensionless tunnel conductance a = Rk/Rt, where Rk = h/e^ is the quantum 
resistance and Rt is the tunnel resistance, predicts that when the thermal energy fc^T is 
much smaller than the charging energy = e^/2C, the conductance of the SET vanishes for 
CgUg/e = n, and the conductance attains a maximum for CgUg/e = {n+ 1/2), where n is an 
integral number When the tunnel resistances Rt are not much smaller than the quan- 
tum resistance, the Coulomb blockade of electron tunneling is softened due to co-tunneling 
and higher order tunneling processes 0. 

In order to describe ferromagnetic single-electron transistors (FSET's), the spin depen- 
dence of the density of states and the transition matrix elements must be taken into ac- 
count. Theoretical works on this subject based on the lowest order sequential tunnelling 
processes have recently been carried out |Tl|-p!^. By including the co-tunneling processes 



in the Coulomb blockade regime, it has been shown [|l4| , p!5| that the TMR is larger in such 
ultrasmall structures than in macroscopic devices, where the charging energy is negligible. 
When the junction resistances are smaller than the quantum resistance, the calculations 
based on sequential tunneling or co-tunneling become inaccurate because higher order tun- 
neling processes determine the transport properties. It is so far unclear if the TMR of the 
FSET is further enhanced by the higher order tunneling processes. This question cannot 
be addressed from the results of the sequential tunneling or co-tunneling formulations. We 
will employ a non-perturbative approach to clarify this issue and show that indeed there 
is a further large enhancement of the TMR in the strong tunneling regime as compared 
with the weak tunneling case where the junction resistances are larger than the quantum 
resistance. This enhancement of the TMR is much larger than what calculations based on 
the co-tunneling processes predict. 

We consider the situation when the two leads are made of the same material and have the 
same magnetization directions, and denote the resistance of the FSET as i?^, where 77 = P 
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(?7 = AP) for the case that the magnetization of the central grain is parallel (antiparallel) 



to the magnetization in the leads [^p!THl^. The magnetoresistance ratio is defined by Eq. 
(|^). In order to calculate the resistance of the FSET for arbitrary tunneling resistances 
the tunneling Hamiltonian must be treated in a nonperturbative manner. This can be 
realized via the path integral approach M. We use the standard single-electron tunneling 



Hamiltonian with spin- dependent density of states and tunneling matrix elements ||11|J14 
and limit our discussions to the linear response regime. 

At sufficiently high temperatures {ksT ^ Ec) the system is in the classical regime, the 
Coulomb charging effects are negligible and the conductance of the system can be read- 
ily found. The dimensionless conductance of the FSET for the parallel magnetic align- 
ments is Op = q;m [1 + (1 — — + + -P')]) where au is the dimensionless 
majority-band conductance of the FSET with parallel magnetic alignments, and the spin 
polarizations of the leads and island are denoted by P and P', respectively. Similarly, 
the dimensionless conductance of the FSET for the anti-parallel magnetic alignments is 
a^p = aM [(1 ~ + -P) + (1 ~ + -P')]- The classical (high-temperature) magne- 

toresistance ratio is 

pel pel o p p/ 

cl _ -ftAP — -ftp ^ zr-r- 
^ ~ Rl,^ 1 - PP' ■ ^ ^ 

At lower temperatures the Coulomb charging effects must be taken into account. In this 
case, the resistance of the FSET can be calculated via Kubo's formula [|T^], 

= in lim u-^Qm{ lim / dre*'^'>(r)} , (3) 

Rjj i^^O iijJi^iO+iS Jo 

where the Matsubara frequencies are uji = 27r/ /jS {I is an integral number) and 



fc=— oo 



oo 

2iTikn„ 



X / Dipe-'^^^'^\{T) cos[ip{T) - ip{0)] . (4) 
Jbk 

where the damping kernel x{^) is an even function with a period f3 and Fourier components 
x{uji) = —\uJi\/4:7i and n^^ = CgVg/e is controlled by the gate voltage. In Eq. the 
boundary condition for the path- integral is 6^ ^ (f{f3) = (f{0) + 2Tik and is the partition 
function of the unbiased device, namely a single-electron box (SEB) Z^, is readily 

represented as a path integral 

oo „ 

Zr,= e^"^'^"- / Dy^e-^"'^! . (5) 

The action is given by 

Sr^[ip\ = dr 

•J 



fir) 



-a^n dr dr'xir - r') cos[v9(r) - (/^(r')] , (6) 
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which depends on the relative orientations of the magnetizations in the central grain and 
the leads through the second term in Eq. (P) describing the contribution of the tunnel 
junctions. If this term is neglected, the charging energy determines the action, and the 
sequential tunneling processes are recovered. If one expands the partition function (|^) in 
terms of the dimensionless conductance, a Taylor series is obtained containing sequential 
tunneling, co-tunneling and higher order contributions. However in practice, it is quite 
tedious to evaluate the path integral term by term in this way. Instead, we will evaluate the 
DC current nonperturbatively. 

Note that the current-current correlation functions of the SET contain both sine and 



cosine terms in the phase differences, as described in Ref. [jTO[. However, using current 
conservation through the first and the second junction it can be seen that the contribution 
from the sine term to the conductance is proportional to the contribution from the cosine 
term and we obtain the simplified expression (^). 

Let us first discuss the most interesting case when the Coulomb charging energy is much 
larger than the thermal energy and there is a strong influence of the Coulomb charging on 
the current. In this case the conductance of the FSET can be varied by tuning the gate 



voltage PJT^l- We use Monte Carlo simulations to compute r(r) (^), and then calculate the 



conductances via the Kubo's formula (^. Since in the Coulomb blockade regime, the TMR 
for integer is at a maximum [^,14|, we will concentrate on this case to find the largest 
TMR and we set ripx = 0. The problem is reduced to the evaluation of expectation values. 



which we will treat by Monte Carlo simulations [|T7 



The computations have been carried out using the standard Metropolis algorithm. The 
convergence of the simulations have been checked by increasing the number of Trotter in- 
dices, which are the numbers of slices in the imaginary time interval [0, /3] and are found 
to scale linear in (3. Therefore at low temperatures where (3 is large the simulations are 
time-consuming. Samplings have been done every 5 passes, and the program has been ran 
for a sufficiently long time to equilibrate the system before sampling. In order to guarantee 
a high precision of the Monte Carlo simulations, the number of samplings have been deter- 
mined so that the result with doubled sampling number 2N is the same as the result with 
sampling number N . The typical number of samplings are found to be of order of 10^ for the 
parameters described below. This procedure is necessary because we calculate the Fourier 
components from the Monte Carlo data and then carry out the analytical continuation 
using Fade approximate to get the conductance of the FSET [p!8ip!9| . 

When the tunnel resistances are close to the quantum resistance, the co-tunneling result 



of the conductance of the SET p,14,20| is recovered. In order to find out if higher order 



tunneling processes will enhance the TMR, we show the result for the conductance of a 
typical value in the strong tunneling case, namely = 10. We consider Ni leads with spin 
polarization P = 0.23, and Co island with P' = 0.35, as in the experiment by Ono et a/ 0. 
For this device, the prediction of the calculations of the TMR in the classical regime (high- 
temperature) yields a value of 17.5 % and in the Coulomb blockade regime the co-tunneling 



formalism gives a value of 38 % |n] . The computed results of the TMR as a function of 



Ec/ksT are shown in Fig. |l[ In the low temperature quantum transport regime, the TMR 
is denoted by 'j^^. We see that the TMR is strongly enhanced by lowering the temperature 
and reaches the co-tunneling result when Ec/ksT ~ 15. In the case of co-tunneling, 
the TMR saturates at this temperature and remains constant at lower temperatures. We 
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would like to emphasize that in our case by including higher order tunneling processes 
a further large enhancement of the TMR occurs at even lower temperatures, e.g. when 
Ec/kBT ~ 40 the TMR is enhanced to 90 %, which is a much larger value than in the 
co-tunneling regime. At even lower temperatures a further enhancement is expected. The 
fact that there is a large enhancement of the TMR in the strong tunneling regime can be 
understood as a consequence of higher order tunneling. In the co-tunneling regime the 
conductance of the device is proportional to the square of the tunnel conductance which 
roughly doubles the TMR with respect to the classical value. In our case with = 10 
the conductance increases even faster than the square of the tunnel conductance and hence 
the TMR is further enhanced. The numerical results of the TMR in the low temperature 
quantum regime can roughly be fitted by 

= + 1^^- (i + ^^1) {EJkBTf , (7) 

where we choose ^"^^ = 0.0006 so that the fitting curve is close to most of the data. Note 
that although //^^ is quite small, the second term on the right hand side of (|^ can be as 
large as unity when the temperatures are sufficiently low leading to a large enhancement of 
the TMR. We have thus shown that there is a large enhancement of the TMR in the strong 
tunneling regime at low temperatures, which is well beyond the results from the sequential 
tunneling and co-tunneling formalisms. 

Let us now consider the semiclassical regime when the thermal energy is larger 
than the Coulomb charging energy Ec- The fluctuation amplitude of the phase variable 
can then be estimated from the kinetic part of the action giving a value of the order of 
[Ec/kBTY^'^, and the dominating part of the action is Gaussian with some corrections from 
the nonquadratic parts of the tunneling term in the total action. By expanding the cosine 
function in the action Eq. (P) and in Eq. (^, the resistance of the FSET to the second order 
in Ec/ksT is 

The calculation leading to the above equation is very similar to the one of the nonmagnetic 



SET carried out in Ref. ||T0[. The two flrst terms in Eq. (H) have also been obtained in 
Ref. [|T^. The third term could also have been derived from the formalism used in Ref. 
14 1 , but the last term is due to higher-order tunneling. It is only this term that gives a 
temperature-dependent contribution to the TMR at high temperatures. The TMR in the 
semiclassical regime is 

^se ^ ^cl ^ ^se ^-^ ^ ^cl^ {EjkBTf , (9) 

where 

/i^^ = 0.04aAf(l + PP')/ [(1 + P){1 + P') ] . 

For P = 0.23 and P' = 0.35, we flnd /i"" ^ 0.013aA/ and the TMR for intermediate 
Coulomb charging energy in the strong tunneling regime, e.g. for E^. = kBT /A and Um = 10, 
increases by less than 1 % giving a TMR much smaller than the low-temperature co-tunneling 
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result of 38 %. The above analysis shows that when the thermal energy is dominating, the 
Coulomb charging effect does not drastically change the TMR, and the role of the higher 
order tunneling processes is not very important. In order to make the TMR significantly 
different from the corresponding values of a macroscopic device, one has to minimize the 
junction capacitances so that the devices work in the low temperature regime. 

It is interesting to note that the enhancement of the TMR the in low temperature 
quantum regime (|^ has a similar temperature-dependence as the enhancement of the TMR 
in the high temperature semiclassical regime (|^), with fi^^ much smaller than /j,^'^. This 
can be understood in terms of smearing of the Coulomb blockade by the strong tunneling 
processes. Taking into account that at low temperatures the renormalized charging energy 
of the FSET, E*, is one order of magnitude smaller than the bare charging energy Ec when 
aM = 10 [ig, we see that {El/ksTf ~ Qm^E^keTf, consistent with Eqs. (0) and (D 
where ~ 0.01. The argument used here will fail at sufficiently low temperatures, 

where the TMR can hardly be described by a power law of {Ec/ksTY . In such an extremely 
low temperature and strong tunneling regime, more investigations need to be done. 

In conclusion, we have derived formulas for the resistance of the ferromagnetic single- 
electron transistor which are valid also in the regime when the junction resistances are 
smaller than the quantum resistance. In the strong tunneling regime we find that the 
magnetoresistance ratio is largely enhanced by the higher order tunneling processes at low 
temperatures. The magnetoresistance ratio in this regime is much larger than the value 
predicted by only including sequential or co-tunneling processes. We find that the TMR 
can at least be enhanced by a factor of 5 relative to the classical value at low temperatures. 
When the thermal energy is larger than the charging energy, the enhancement of the TMR 
is small, even when the junction resistances are much smaller than the quantum resistance. 

The authors are grateful to G. E. W. Bauer, K.-A. Chao, J. Inou, A. N. Korotkov, 
Yu. V. Nazarov and A. A. Odintsov for stimulating discussions and R. Egger for valuable 
support of the Monte Carlo simulations. X. H. W. would like to thank the Swedish TFR 
for financial support and A. B. acknowleges financial support from the EU TMR Research 
Network No. FMRX-CT96-0089 (DG12-MIHT). 



6 



REFERENCES 



^ also at Philips Research Laboratories, Prof. Holstlaan 4, 5656 AA Eindhoven, The 
Netherlands. 
[1] M. JuUiere, Phys. Lett. 54 A, 225 (1975). 

[2] J. S. Moodera, L. R. Kinder, T. M. Wong and R. Meservey, Phys. Rev. Lett. 74, 3273 
(1995); J. S. Moodera and L. R. Kinder, J. Appl. Phys. 79, 4724 (1996). 

[3] P. M. Levy and S. F. Zhang, Current Opinion in Solid State & Materials Science 4, 223 
(1999). 

[4] K. One, H. Shimada and Y. J. Ootuka, J. Phys. See. Jpn. 66, 1261 (1997). 

[5] L. F. Schelp, A. Pert, F. Fettar, P. Holody, S. F. Lee, J. L. Maurice, F. Petroff and A. 
Vaures, Phys. Rev. B 56, R5747 (1997). 

[6] S. Sankar, B. Dieny and A. E. Berkowitz, J. Appl. Phys. 81, 5512 (1997). 

[7] D. V. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 65, 2446 (1990). 

[8] D. V. Averin and K. K. Likharev in Mesoscopic Phenomena in Solids, edited by B. 
L. Altshuler, P. A. Lee, and R. A. Webb (Elsevier, Amsterdam, 1991); G.-L. Ingold 
and Yu. V. Nazarov in Single Charge Tunneling, Vol. 294 of NATO Advanced Study 
Institute Series B, edited by H. Grabert and M. H. Devoret (Plenum, New York, 1992). 

[9] G. Schon, and A. D. Zaikin, Phys. Rep. 198, 237 (1990) and references therein. 
[10] X. H. Wang, Phys. Rev. B 55, 12 868 (1997). 
[11] J. Barnas and A. Pert, Phys. Rev. Lett. 80, 1058 (1998). 

[12] A. Brataas, Yu. V. Nazarov, J. Inoue and G. E. W. Bauer, Phys. Rev. B 59, 93 (1999); 

Eur. Phys. J. B 9, 421 (1999). 
[13] A. Korotkov, Phys. Rev. B 59, 89 (1999); Superlatt. and Microstr. 25, 259 (1999). 
[14] S. Takahashi and S. Maekawa, Phys. Rev. Lett. 80, 1758 (1998). 
[15] H. Imamura, S. Takahashi and S. Maekawa, Phys. Rev. B 59, 6017 (1999). 
[16] M. Biittiker, Phys. Rev. B 36, 3548 (1987). 

[17] X. H. Wang, R. Egger and H. Grabert, Europhys. Lett. 38, 545 (1997); and references 
therein. 

[18] G. D. Mahan, Many-Particle Physics, Chapter 3, (Plenum Press, New York, 1990). 
[19] S. Chakravarty and J. Rudnick, Phys. Rev. Lett. 75, 501 (1995); and J. E. Hirsh, in 

Quantum Monte Carlo Methods in Equilibrium and Nonequilibrium Systems, edited by 

M. Suzuki (Springer, Berlin Heidelberg, 1987). 
[20] P. Joyez, V. Bouchiat, D. Esteve, C. Urbina and M. H. Devoret, Phys. Rev. Lett. 79, 

1349 (1997). 



7 



FIGURES 

FIG. 1. The tunnel magnetoresistance ratio (TMR) as a function of the dimensionless charging 
energy in the strong tunnehng regime for aM = 10. The diamonds arc Monte Carlo data with 
errors smaller than the sizes, and the solid line is calculated via Eq. 7 for fi"^" = 0.0006. 
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